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ABSTRACT 

In the entrance region of a nonequilibrium MHD generator one finds an abrupt rise 
in electron temperature followed by a four- or five-decade increase in electron density. 
We consider the stability of this region of abrupt electron density rise. A solution is 
found in the limit that the ambipolar diffusion field is small compared to the transverse 
applied electric field. The solution is found to be stable in the limit of zero magnetic 
field but to be unstable in the presence of a magnetic field. 



ELECTROTHERMAL INSTABILITIES IN THE ENTRANCE REGION 
OF AN MHD GENERATOR 
by J. Marlin Smith 
Lewis Research Center 

SUMMARY 

Previous analyses of electrothermal instabilities in MHD generators have considered 
fluctuations from a steady state solution pertinent to the infinite homogeneous medium. 

In the entrance region of a nonequilibrium generator this condition is not valid since in 
this region one finds an abrupt rise in electron temperature followed by a four- or five- 
decade increase in electron density. We consider the stability of the subregion of abrupt 
electron density rise. The steady state solution is one of uniform gas dynamic proper- 
ties and electron temperature but with an electron density gradient which is solely a 
function of the ionization rate. This rate is evaluated at an electron temperature ob- 
tained by equating the Joule heating to the elastic collisional and ionization energy loss. 
The governing equations upon which the perturbation is taken are those of reference 1 
with the exception that recombination is neglected in the electron continuity equation. A 
solution is found in the limit that the ambipolar diffusion field is small compared to the 
transverse applied electric field. The solution is found to be stable in the limit of zero 
magnetic field but to be unstable in the presence of a magnetic field. 


INTRODUCTION 

In contrast to experiments performed with equilibrium MHD generators, nonequilib- 
rium generators have conspicuously failed to operate at anything near the performance 
predicted by theory. These negative results have been attributed to such factors as Hall 
current leakages resulting from the higher Hall parameters required to operate in the 
nonequilibrium mode, insufficient residence time of the working fluid in the generator 
for the nonequilibrium ionization to buildup to its steady state value, instabilities arising 
from the nonequilibrium plasma state, etc. In this report another possible source of 
difficulty is considered. This is the possibility of the formation of instabilities in the 



ionization region of the entrance of the generator which would inhibit the growth of ion- 
ization to the fully nonequilibrium level which present theory now assumes. 

Previous analyses of electrothermal instabilities in MHD generators have been con- 
cerned with fluctuations occurring in regions well removed from any large gradients of 
gas dynamic or electrical properties. In this regime the steady state solution is that 
pertinent to the infinite homogeneous medium or small deviations therefrom. In the en- 
trance region of a nonequilibrium generator these solutions are only valid for extremely 
short wavelength component in the gradient direction since in this region one finds an 
abrupt rise in electron temperature followed by a four- or five-decade increase in elec- 
tron density occurring within a few centimeters. It is therefore the stability of fluctua- 
tions with wavelength component the order of the ionization region which is of concern in 
this report. 

The analysis is limited to the segmented electrode generator configuration operating 
in the Faraday mode. The ionizing electric field is restricted to that induced by the gas 
flow through the magnetic field or this field supplemented by an electric field applied 
solely in the Faraday current direction. All equations are in SI units and all symbols 
are defined in the appendix. 


STATEMENT OF PROBLEM 
Model 

To facilitate the analysis a simplified model of the ionization region based upon the 
numerical results of Bertolini (ref. 2), is taken. In the model considered, it is assumed 
that at the entering electron temperature (gas temperature) the ionization rate is small 
compared to the electron heating rate so that the entrance region can be considered to 
consist of two subregions. In the first subregion the electron temperature is abruptly 
elevated with negligible change in electron density until the point is reached at which the 
Ohmic heating is balanced by elastic electron- neutral atom collisional and inelastic ion- 
izing collisional energy losses. This is followed by a second subregion in which the 
electron density rapidly rises at constant electron temperature. 

This model is realistic due to the very strong exponential dependence of the ioniza- 
tion rate on electron temperature so that the distance over which the electron density 
builds up and reaches equilibrium is controlled by the peak electron temperature. That 
the electron temperature remains constant in the second region is a consequence of the 
electron density dependence of the electron heating equation and will be demonstrated at 
a later point in the analysis. 
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Assumptions 


On the basis of the previous model this report is limited to the study of the stability 
of fluctuations occurring in the second subregion. Furthermore, as in the case of homo- 
geneous electrothermal instabilities, the following limitations on the analysis are made: 

(1) Throughout the analysis the gas dynamic properties, that is, gas density, veloc- 
ity, and temperature, are taken to be constant. In the steady state solution this is a 
valid assumption since these properties change on a scale which is the order of the gen- 
erator length (typically the order of meters) while the ionization region is ideally the 
order of centimeters. Fluctuations in the gas dynamic properties are ignored relative 
to those of electron density and electron temperature due to the fact that the relatively 
lighter and more mobile electrons respond to disturbances of greatly different frequen- 
cies than do the heavier and less mobile neutral gas atoms which dominate the gas prop- 
erties. Furthermore, due to the low degree of ionization in the region of interest, the 
gas dynamic properties are weakly coupled to the electrical properties. 

(2) Only the propagation of magnetohydrodynamic type disturbances are considered. 
This restriction embodies the following limitations: 

(a) The displacement current is neglected in Maxwell's equations. Under this 
restriction the relatively high frequency transverse electromagnetic disturbances are 
not considered. 

(b) Induced magnetic fields are ignored, that is, the analysis is limited to con- 
ditions of low magnetic Reynolds number which is, in general, satisfied in MHD gener- 
ator applications. 

(c) The plasma is assumed to be quasi- charge neutral, that is, to zero order 
the electron and ion number densities are equal. 

(3) Terms of the order of the electron to heavy particle mass ratio are neglected. 

(4) Ion slip is neglected. This is justified since in cases of interest the Hall param- 
eter is less than 10. 

(5) Consider only propagation in the plane perpendicular to the applied magnetic 
field. 


Equations 

For the previous model and assumptions the analysis is completely determined in 
terms of the electron density and electron temperature by the following set of equations. 

Maxwell equations . - Under the previously discussed conditions of charge neutrality 
and neglect of the induced magnetic field these equations reduce to 
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( 1 ) 


V X E = 0 

V.j = 0 (2) 

Electron continuity equation . - In the region of interest the dominant process is that 
of ionization and the analysis is restricted to the case in which the ionization process is 
dominated by electron-neutral atom ionizing collisions. Then 


_a 

0t 


n 


+ v 




V "e = V Vi 


(3) 


It is noted that the equation depends upon the gas flow velocity Vq rather than the elec- 
tron velocity v . This arises as a result of the neglect of ion slip and the assumption of 

V ■ 

quasi- charge neutrality for which the definition of J_ is n g e(vQ - v ) and V • j. = 0. 
Therefore 


V* 


n,v 

e-e 


V- 




"e^O = —0 


Vn, 


(4) 


that is, gas flow velocity is constant. 
Generalized Ohm’s law. - 


j = <xE* - — j x B n + — VP (5) 

"VO "Vo 

where E* = E + Vg x Bq is the electric field in the frame of reference moving with the 
velocity v Q and Bq is the applied magnetic field. 

In the above form of the generalized Ohm’s law the contribution from the heat con- 
duction is neglected. Its effect is to contribute terms of the form Vn and VT which 

_ v C 

lead to a slight modification of the VP term already appearing in the Ohm’s law equa- 

V 

tion. However, the slightness of the correction does not seem to justify the additional 
complexity resulting from its inclusion and therefore is ignored. 

The electrical conductivity a is defined in the usual manner as 


a = 


n e e 


( 6 ) 


"Vo 


where, as in reference 1, the total electron momentum collision frequency is taken 
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to be constant throughout the analysis. This condition is justified, particularly in the 
region considered here, because the electron density is small so that the collision fre- 
quency is dominated by electron- neutral atom collisions. In this case the frequency de- 
pends upon the neutral atom density (=* gas density) which by prior restriction is constant. 
We need then only require that the energy dependence of the cross section be such that 
the frequency have a negligibly small dependence upon electron temperature (nearly 
Maxwellian molecule interaction). This is assumed to be the case in this analysis. 

Electron energy equation . - 



+ n el0 




= " 2 ir 1 '0 n ef k(T e - T 0> + 1 ' 5* " ( !cT i +| kT e)(|.“e + ^0 ’ Vn e) ™ 

where the VT g contribution to the heat conduction term and radiation losses are ne- 
glected. These terms contribute wavelength dependent damping terms to the dispersion 
relation which decrease with increasing wavelength relative to the wavelength indepen- 
dent elastic collision damping. The analysis is therefore restricted to wavelengths of 
sufficient magnitude so as to make these terms ignorable. For a further discussion of 
this point see reference 3. 


ANALYSIS 


The analysis is carried out for a segmented electrode MHD generator operating in 
the Faraday mode. The coordinates are oriented such that the x-axis is in the fluid flow 
direction and the z-axis is in the direction of the applied magnetic field. 


Steady State Solution 

Under the previously discussed assumptions the zeroth order or steady state solu- 
tion is depicted by a region of constant fluid dynamic quantities and a constant electron 
temperature which is given by equating the Joule heating to the elastic and ionization 
collisional energy losses. From equation (7) the steady state electron temperature is 
then given by 
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( 8 ) 


HI 


io • lo = 2 ^ VeO ~ k < T eO " T 0> + ( M i + § kT eO )Zo ' Vn , 


eO 


where the subscript ”0” refers to the zeroth order or steady state solution. 
From equation (3) the steady state electron density is given by 


lo * Vn e0 _ n eO n aO l/ iO 


(9) 


Since the fluid properties are constant n & Q ^ constant and since only depends upon 
the constant electron temperature it is constant. Therefore, with Vq a constant and 
orientated along the x-axis 


v n • Vn„ n = v n — In n eQ = n a0 ^ iQ = constant 


n 


eO 


^0 


eO 


dx 


( 10 ) 


Obviously the space dependence of the steady state electron density is exponential with 
an e-folding length X. defined by 


1 d WiO 

— = — In n 0 = ■■ 

h * v 0 


(ID 


Under conditions typical of nonequilibrium generator operation the incoming electron 
density which is in Saha equilibrium with the electron temperature equal to the gas tern- 

1 fi 

perature is of order 10 electrons per cubic meter. In the ionization region the elec- 
tron temperature is elevated by the preferential Joule heating of the electrons to a level 

at which the electron density in Saha equilibrium with the electron temperature is in the 
20 21 

range of 10 to 10 electrons per cubic meter. Furthermore, if this increased ioniza- 
tion is to be effectively used in the generator then the ionization region must be the order 
of a few centimeters at most. Therefore, the e-folding length must be on the order of a 
tenth of a centimeter. 

The steady state Ohm's law relation as obtained from equation (5) is 

e . kT P o - 

lo = a 0^0 " ~~"!o x ®0 + a 0 ~~~ V ln n e0 ( 12 ) 

m e v 0 e 

With the definition of the segmented Faraday generator, 
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3 0x " 0 


( 13 ) 


and with the usual insulating walls in the magnetic field direction so that 

j 0z - E Sz = 0 < 14 > 

the other components of equation (12) are 


where 


E 0x E 0x ^0^ E 0y " V 0 B 0> 


kT 


eO d 


In n 


e dx 


eO 


j 0y = CT 0 (E 0y “ V 0 B 0 ) 


*0 



in 

e 0 


(15) 

(16) 

(17) 


is the Hall parameter. 

Equations (8), (10), and (13) to (17) give a complete description of the steady state 
solution in the region of interest. Before proceeding to answer the question as to 
whether or not this solution is stable to small perturbations we will show that the as- 
sumption of constant electron temperature is consistent with this set of equations. 

Since there is no charge separation and Vq and Bq are constant 

E 0y = E 0y _ V 0 B 0 = constant < 18 ) 

Then dividing equation (8) by <7 q and using equation (16) 


j HI 

a 


Y f k < T e0 - T 0) + ( kT i +1 2 kT eo) IQ ' ‘ V ln n e0 


(19) 


Since by prior argument Vq and Tq are constant, it is obvious from equations (11) 
and (18) that equation (19) requires T g Q to be constant. 

It is also convenient to note from equations (10), (17), and (18) that equation (15) 
requires that 


Eqx = constant 


( 20 ) 
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and, therefore, that the steady state electric field 


Eq = constant 


( 21 ) 


Furthermore, since 


<r 0 = 


n 


m e v 0 


eO 


the conductivity has the same exponential x-direction dependence as n g g so that by 
equations (18) and (16) jgy must have this same dependence. The steady state solution 
can now be summarized as 


n e0 ~ e 


(iAj)x 


Eq x , Eg y = constant 


E 3z = 0 


■'Ox ■'Oz ® 

► 

(lApx 
j 0y~ e J 


T e0 = constant 


( 22 ) 


(23) 


(24) 


(25) 


Perturbation Equations 

The essence of this analysis is to determine whether an infinitesimal perturbation 
of the static system will initially grow or decay. To this end we therefore consider 
solutions of the governing equations of the form 

f = f 0 + f’(r,t) (26) 

where fg is the steady state solution and f’ is an infinitesimal perturbation to this 
solution such that products of the f’ can be neglected. Substituting solutions of this 
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form into the governing equations, neglecting perturbations in the gas dynamic quantities, 
and utilizing the steady state solutions the equations for the perturbation quantities be- 
come 

Maxwell’s equations 


Electron density 


VXE' = 0 

v. r = o 



at n 


eO 




VI 


v i0 


lo • V ln n e0 


(27) 

(28) 


(29) 


v • is determined as a function of T by approximating the ionization cross section by a 

X v 

linear curve fit to the low energy portion of the curve and averaging over a Maxwellian 
electron distribution so that 


V . ■*** 

1 


. 1 / 2 , 


1 +2 


- T i/ T t 


(30) 


where Tj is the temperature below which the ionization cross section goes to zero. At 
low electron densities where multistage ionization processes are negligible, T. is the 
temperature corresponding to the ionization potential. This is assumed to be approxi- 
mately true in the region being considered here. The perturbed value of equation (30) is 
then 


v 


t 

i 


r^> 


V 


iO 



(31) 


where we have neglected terms the order of T 6 q/T j since in the regime of interest 
T eO /T i « 1 . Equation (30) then becomes 


n' 


9t n 


+ v, 


eO 



(v Q • V ln n e0 ) 


Tl 


eO 


(32) 
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Generalized Ohm's law 



Electron energy equation 



Solution of the Perturbation Equations 

We first write equation (28) in the form 



It is now noted that equations (27) and (32) to (35) are simply a linear set of coupled 
equations with constant coefficients for the variables 



Therefore, the solution can be taken in the form 
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rs 


ro 






n. 


n 


eO 


a 0 


i, 



e i(Z-r+wt) 


(36) 


iv yv/ r>j 

where E, n e / n e Q> (1 /oq) j , and T e /T g Q are all constant and considering only propaga- 
tion in the plane perpendicular to the applied magnetic field (Z , Z , Z ) = (Z , Z„, 0). 

a y .4 a y 

It should be noted at this point that the normalization of the perturbation equations 
has resulted in solutions in the form of superimposed waves. However, this differs 
from the homogeneous ease in that not all of the perturbed quantities have a wave- like 
form. Consider, for example, the perturbation in the electron density which from equa- 
tions (22) and (36) is of the form 


n' 


(lAi) x i(£. r+cot) 


This solution exhibits growth regardless of the eventual solution of the dispersion rela- 
tion. However, it is noted that this growth is at the same rate as the zeroth order elec- 
tron density so that if the perturbation is small to begin with it will be small to the same 
ratio throughout the region. This is not unstable behavior since the infinitesimal per- 
turbation remains infinitesimal relative to the zeroth order. Therefore, just as in the 
case of the infinite homogeneous medium, when the imaginary part of w is negative 
the perturbation grows without bound and the system is termed unstable; a negative value 
of the imaginary part of w in the case considered here indicates that the perturbation 
will grow faster than the zeroth order and hence that the system is unstable. 

Introducing the solutions represented in equation (36) into equations (27) and (32) to 
(35) we obtain 


1x1 = 0 


(37) 


o>* — — 
n e0 



(X0* V ln n e0> 




<->*/ 



+ E - 


e 


"Vo 



. kT e0 

x Bp + i— 



(38) 


+ 


e \ n e0 



(39) 
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where co* = co +_Z • Vq is the frequency in the frame of reference moving with the gas. 

The dispersion relation can then be determined from the condition for the existence 
of a solution, namely, that the determinant of coefficients be zero. However, it is 
somewhat more illustrative to derive the relation by the process of eliminating the 
various variables in the electron energy equation (eq. (40)). Therefore, the procedure 
to be used in this report is to use equations (37), (39), and (41) to eliminate the 



term in terms of T e /T eQ and n /n eQ , then eliminate n e /n eQ by equation (38). 

The simplest way to do this is to observe from equation (5) that j • E* can be 
written as 




1 

n e e 



= — - j- j - — iV vp 

a 0 a ~~ n e e \ a 0 7 

which upon linearization and substitution of equations (12) and (36) yields 
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rV /Ny-I 

To determine (1/ctq) j, the cross product |xux (1/ctq) j J is formed from equation (39), 
which, with the use of equations (37) and (41) yields 



(43) 


The quantity 


— j I • V In n 


i o 


0 


eO 


is then determined from equation (43) by dotting V In n q into the equation to obtain 



rir vlnn e 0 = 

i Or 


i4 + e x - (3 oVH 


I 2 +i(! x - (3 0 ! y ) 1 


X i 


n^ T 0 /kT Qn dlnm n 

n e0 i e0 


dx 


(44) 
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Substitution of equations (43) and (44) into equation (42) after considerable algebraic ma 
nipulation yields 


1 ■!)•■;+ 



rsj 

E s 



F * n e . H e0 

— o 1_ — 

n e0 e 



eO 


a - iA 


n. 


g - iG 


4 . <*x - W/ ne0 ,4 _ ‘ e0 


,2 T, 


l + 


r + 




(45) 


where using equations (15) and (16) to eliminate Eq x and i 0y in terms of E^. 


J 0y 


a - 2 l 2 l y (l y + P 0 l x )E $ - l\(l x - Vy>E5y(-^ ^ 0l 


(46) 


rkT ~ d In n 


A o 2 iy (i y + P 0 l x )(,l K - e 0 9 A E$ + y 4 A Ej y ( 40 e0 ' 


e dx 


(47) 


g = 2 




El I ^eO i^eo\ + j2 A^eO H^eO j (48) 

i dx / M e dx 


J 0y 


« < /kT on d In n. n \ 0 , /k'T n dlnn 

G - 2i y (i y + P 0 l x ) A Eji-£2 -^) + l]d x - Vy> A | — — — ^ (49) 


Xjl e dx 


Upon substitution of equation (45) into equation (40) and using equation (8) to eliminate 
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we obtain 
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which is the general dispersion relation for the problem under consideration. Unfortu- 
nately it is too complicated to consider for other than numerical evaluation. Fortunately, 
however, in the regime of parameters of interest in most nonequilibrium MHD experi- 
ments the relation can be greatly simplified. 

It is noted that the coefficients a, A, g, and G are functions of the Faraday elec- 
tric field Epy and the ambipolar electric field 

kT eO d ln n e0 
e dx 

and by virtue of equations (16) and (19) so are the coefficients of equation (51). Under 

O 

typical generator operating conditions E^ =* 2x10 volts per meter and T 0 q ^ 3000 K. 

2 3 

At this temperature the ambipolar field ^3x10 volts per meter (d In n e0 /dx-10° 1/M, 
see discussion following eq. (11)). Therefore 

^eO dIn n e0 
_£ — « i 

TT1 * 

E 0y 

The coefficients of equation (51) are now simplified by retaining only the dominant power 
of the above ratio in each of the coefficients to yield 


w 


*2 


i + i ^L 


3 ,4, V V/ 
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I* in . m e 

- — — + i s 2 Vq + 

n e0 e m a 


1 + — /— + 


3 \ 2 T eo /\ 2 T 




eoy 


Zb- 


V In n 


eO i 
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= 2 






+ l- 




z 4 + (Z x ~ 


A “ ^0 V 

L H — — — — — 


4 


2 V + - (- + — L] v . v In n 0 

. m a 0 n* T J 


<! T.\ _ 

_ + Vo-Vlnn^ (52) 


[2 T 


e0/ 


16 



10 — 
It is further noted that typically Vq =* 10 collisions per second and m e /m a ^ 10“ 

so that 


111 c i 

2 i/ Q « 10° — (53) 

m_ see 

ci 


Also 


and 


so that 


v n 10 3 JL 
u sec 
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d ln n eQ 
dx 



(54) 


Furthermore 


iio 


n e0 e m e^0 


*y E 0y — ^0^y v 0 ~ % 


2ir 


where Xy is the wavelength in the y-direction. In general, < 1 so that for wave- 
lengths the order of centimeters in the y-direction 


l ‘ jn c 

10 & (55) 

n e0 e 


From relations (53) to (55) we conclude that the factor 


d In n. 



17 



is by far the dominant one in equation (52) and hence we seek a solution in the limit as 
this factor approaches infinity. For this limit, the real and imaginary parts of co be- 
come 




* 

w r = 


,4 , ~ V/ 


- v 0 -Vlnn e0 


(56) 


*2 

X i 


lljl. v+W 

W i = 2 — 7 V ^ ne0 

,4 + (jx ~ W 


(57) 


X 2 

X i 


That this is the correct solution to equation (52) to the order of the ratio of relations 
(53) and (55) to relation (54) can be shown by direct substitution of the solution back into 
equation (52). 

It is noted from equation (57) that for zero magnetic field, |3 q = 0, ok is positive 
and hence the system is stable at least up to terms the order of the ratio of ambipolar 
to Faraday electric fields. However, with /3 q > 0, ok becomes negative for orienta- 
of the wave vector for which 


0 0 — <-l (58) 

V 

Since all orientations are possible, the result indicates that the system is unstable for 
all nonzero values of magnetic field. However, again the result is only correct up to 
terms the order of the ratio of ambipolar to Faraday electric fields. If these higher 
order terms are positive then there will be a critical value of |3 q for which the system 
goes unstable. 


CONCLUSIONS 

It is concluded that the ionization region of a nonequilibrium, segmented electrode 
MHD generator operating in the Faraday mode is unstable in the presence of a magnetic 
field. To the present order of solution, which is the dominant order in the ratio of the 
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ambipolar to Faraday electric fields, the system is unstable to any nonzero magnetic 
field. It therefore appears that some auxiliary means of ionizing the gas in this region 
may be required or the region should be shielded from the magnetic field. Under in- 
vestigation is the possibility of ionizing by applied electric fields within the magnetic 
field region but in the magnetic field direction or in the Hall field direction. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 7, 1968, 

129-02-08-05-22. 
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APPENDIX - SYMBOLS 


A defined by eq. (47) 

a defined by eq. (46) 

magnetic field vector 
electric field vector 
e unit electric charge 

G defined by eq. (49) 

g defined by eq. (48) 

j electric current density 

k Boltzmann constant 

l wave vector 

m particle mass 

n particle number density 

P pressure 

T temperature 

v velocity 

/3q Hall parameter 


X wavelength or e-folding length 

v collision frequency 

a electrical conductivity 

a) wave frequency 

Subscripts: 

a gas atoms 

e electrons 

i ionization value 

x, y, z components of x, y, z coordinate 

system 

Superscripts: 

* evaluated in frame of reference 

moving with gas 

* perturbed quantities 
defined by eq. (36) 
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